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Abstract 

We complete the recalculation of the available two-loop expressions for the reaction 
77 — > TTTT in the framework of chiral perturbation theory. Here, we present the results 
for charged pions. The cross section and the values of the dipole polarizabilities agree 
very well with the earlier calculation, provided the same set of low-energy constants 
(LECs) is used. With updated values for the LECs at order j?'^, we find for the dipole 
polarizabilities (ai — /3i)7r± = (5-7 it 1.0) x 10~^fm^, which is in conflict with the 
experimental result recently reported by the MAMI Collaboration. 

PACS: ll.30.Rd; 12.38.Aw; 12.39.Fe; 13.60.Fz 
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1 Introduction 



We evaluate the amplitude for 77 'k^'k~ in the framework of chiral pertur- 
bation theory (ChPT ) [1,2,3] at two-loop order, and compare the result with 
the only previous calculation performed at this accuracy [4]. We employ the 
calculational techniques outlined in our previous work [5], which allow us to 
provide a rather compact and easy to use integral representation for the full 
amplitude. We find that the cross section for the reaction 77 — > tt+vt" and 
the dipole polarizabilities agree very well with the results reported in Ref. [4], 
provided that the same set of LECs is used. With updated LECs at order 
[6,7], we find for the dipole polarizabilities the value 
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(5.7 ± 1.0) X 10-^fm^ 



(1.1) 



The MAMI Collaboration [8] has recently reported the experimental result 



The index "mod" denotes the uncertainty generated by the theoretical mod- 
els used to analyse the data. The ChPT calculation is in conflict with this 
prediction, see also [9] for a recent discussion. 

There are good reasons to believe that the chiral prediction Eq. (1.1) is rather 
stable against contributions from still higher orders in the chiral expansion. 
The main point is that the Born-term subtracted amplitude - in terms of 
which the polarizabilities are defined - is regular at the Compton threshold. 
Therefore, the influence of chiral logarithms is suppressed as compared to 
the amplitude at, say, the threshold for nn production, where the amplitude 
generates a branch point. Further, we have analysed the potential influence of 
resonance exchange and could not identify any significant contribution from 
nearby resonances to the combination (cti — f3i)^±. We conclude that the low- 
energy constants at order arc expected to play a negligible effect here. In 
view of these observations, the discrepancy between the chiral prediction and 
the recent MAMI data cannot be explained. 

The article is organized as follows. In Section 2 we spell out the kinematics 
of the process 77 —>■ tt+tt". To make the article selfcontained, we summarize 
in Section 3 the necessary ingredients of the effective Lagrangian framework. 
Here, we also fix the LECs to be used in numerical calculations. In Section 4, we 
display the Feynman diagrams and discuss shortly their evaluation. Section 5 
contains a concise representation of the two Lorentz invariant amplitudes that 
describe the scattering matrix element. Section 6 contains explicit expressions 
for the dipole and quadrupole polarizabilities valid at next-to-next-to-leading 
order in the chiral expansion, together with a detailed numerical analysis and a 
comparison with the recent MAMI data [8] , and with an evaluation from data 
on 77 Ti^Ti~ [10]. The summary and an outlook are given in Section 7. A 
detailed comparison with the earlier calculation [4] is provided in Appendix A, 
and several technical aspects of the calculation are relegated to additional 
Appendices B-D. 



2 Kinematics 

The amplitude describing the process 77 — > ti^ti~ may be extracted from the 
matrix element 



(11.6 ± l.Sstat ± 3.0syst ± O.S^od) X 10-^ fm^ . 



(1.2) 
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( 7r+(pi)7r-(p2) out I 7(?i)7(?2) in > = i{27r)'S\Pf - P,) T^t--^-", (2.1) 
where 

W^;r"^"" =^ / cixe-^(^i^+^^^) ( 7r+(pi)7r-(p2) out | Tj,{x)jM I ) . (2.2) 

Here denotes the electromagnetic current and a = e'^/Air ~ 1/137 is the 
electromagnetic coupling. It is convenient to change the pion coordinates ac- 
cording to (tt^jTT*^) — > (tt^, TT^, TT^) and instead of 7r''"7r~-production, we con- 
sider in the following the process 77 — > tt^tt^, with 

w^jr"^'^" = -w^;r"'"' = -^^^ , (2.3) 

where the relative minus sign stems from the Condon-Shortly phase conven- 
tion. [We use the same sign convention as Ref . [4] .] The decomposition of the 
correlator V^i, into Lorentz invariant amplitudes reads 

V;,^ = A(s, t, u)Ti^^ + B(s, t, u)T2^,, + C(s, t, u)T^^^ + D{s, t, u)Ti^^ , 
Tifj,i, — -s g^j, — qii,q2^ , 

T2nu = 2 s A^A^ -v^ g^^-2v {qiu^^ - q2n^v) , 
Ts/xu = qiij,q2v , 

Tii^u = s{qif,A^ - g2i/A^) -v{qif,qiu + q2^,q2v) , 
A^ = (Pi-P2)m, (2.4) 

where 



s = (?i + ?2)^ ^ = (Pi - qif, H = (P2 - qif, u^t-u (2.5) 

are the standard Mandelstam variables. The tensor V^^, satisfies the Ward 
identities 

q^V^, = q'2V^,.^0. (2.6) 



The amplitudes A and B are analytic functions of the variables s, t and u, 
symmetric under crossing {t,u) — > {u,t). The amplitudes C and D do not 
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contribute to the process considered here, because Si- Qi — 0. 
It is useful to introduce in addition the hehcity amphtudes 



H++^A + 2{iM^-s)B, H^_^^AML_^b. (2.7) 

s 

The hehcity components and correspond to photon hehcity dif- 

ferences A = 0,2, respectively. With our normalization of states (pi|p2) = 
2 (27r)=^p° 5(^)(pi - P2), the differential cross section for unpolarized photons 
in the centre-of-mass system is 



7 -"yy — >'7r~'~7r 2 

-^(^i')H{s,t), H{s,t)^\H^^\' + \H^^\\ (2.8) 



with /3{s) = — AMI/ s. The relation between the hehcity amplitudes M+-|- 
in Ref. [10] and the amplitudes used here is 

M++(s, t) = 2TTaH++{s, t) , M+_(s, t) = 167rQ;S(s, t) . (2.9) 

In the centre-of-mass system, qi-\-q2 — 0, one has qi-pi — \qi\\pi \ cos 9, where 
9 is the scattering angle. Then the Mandelstam variables are given by 

s^A\q\'', t^Ml-{s/2){l- (3{s)cos9). (2.10) 

For comparison with experimental data, it is convenient to present also the 
total cross section for the case having | cos6'| less than some fixed value Z, 



2 * 
an 



(t(s; I cos < Z) = / dtH{s, t) (2.11) 
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t- 



with t± = Ml- (s/2) (1 T (3{s)Z). 



3 The effective Lagrangian and its low-energy constants 



The effective Lagrangian consists of a string of terms. Here, we consider QCD 
with two flavours, in the isospin symmetry limit m„ = = rh. At next-to- 
next-to-leading order (NNLO), one has [2] 
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Aff = ^2 + A + A . (3.1) 

The subscripts refer to the chiral order. The expression for is 



_p2 

D^U = 9^C/ - iiQU - UQ)A^ , Q = |diag(l, -1) , (3.2) 

where e is the electric charge, and denotes the electromagnetic field. The 
quantity F denotes the pion decay constant in the chiral limit, and is 
the leading term in the quark mass expansion of the pion (mass)^, ~ 
M^(l + 0(m)). Further, the brackets (. . .) denote a trace in flavour space. In 
Eq. (3.2), we have retained only the terms relevant for the present application, 
i.e., we have dropped additional external fields. We choose the unitary 2x2 
matrix U in the form 



TT 



/ ^0 



U^a + in/F , a + — = 12x2, tt 



^V2 7r- -7r° 



(3.3) 



The Lagrangian at NLO has the structure [2] 



7 3 _ 7 

A = E li^i + E hiKi = J {D^U D^U^'f + ■■■, (3.4) 

i=l i=l ^ 

where li,hi denote low-energy couplings, not fixed by chiral symmetry. At 
NNLO, one has [11,12,13] 



57 

i=l 

For the exphcit expressions of the polynomials K^, and P^, we refer the 

reader to Refs. [2,11,12,13]. The vertices relevant for 77 — > 7r"'"7r~ involve 
Zi, . . . , from £4 and several q's from Cq, see below. 

The couplings U and q absorb the divergences at order and p^, respectively. 



k^{^^cY-'{ll{^Ji,d) + ^,k} , 



c,: = 



A: 



{cK/., d) - 7f ) A^ - + 7f ^ (/., ^)) A} , 
^^^^,lnc=-l{ln4. + r'(l) + l}. 



(3.6) 
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The physical couphngs are 4) and c[(/i, 4), denoted by in the fol- 
lowing. The coefficients 7.4 are given in [2], and 7^^'^'^"* are tabulated in [12]. 
In order to compare the present calculation with the result of [4], we shall use 
the scale independent quantities k introduced in [2], 

II = 3^ ih + I) , (3.7) 
where the chiral logarithm is / = ln(M^//x^). We shall use [6] 

= -0.4 ± 0.6 , [2 = 4.3 ± 0.1 , [3 = 2.9 ± 2.4 , U = 4:A ± 0.2 , (3.8) 

and 

= - [5 = 3.0 ± 0.3 (3.9) 
obtained from radiative pion decay to two loop accuracy [7,14]. 
The constants occur in the combinations 



al = -40967r^ (6 + c^g - c^^ - 3 0^4 + C35 + 2 c^g - 4 0^7 + C50) , 
= 2567r^ (8 c^g - 8 c^o + C31 + c^2 - 2 0^3 + 4 0^4 + 8 c^q - 4 c^J , 

6'- = -1287r^(c^i + c^2-2c^3-4c^4) ■ 

Their values have been estimated by resonance exchange e.g. in Ref. [4]. We 
have repeated that analysis. Taking into account p, ai and bi exchange which 
contribute with a definite sign, we obtain 

[a\, al, b") = [- 3.2, 0.7, 0.4) [ present work ] . (3.10) 

Unless stated otherwise, we will use these estimates at the scale 11 = Mp. Con- 
tributions from scalar and tensor exchange are of a similar order of magnitude 
(see Table 2 in Ref. [15]). On the other hand, the sign of these contribution is 
not fixed. In Ref. [16], a large Nc framework and the ENJL model were used 
to pin down these constants, with the result 



(a\, al, = 8.7, 5.9, 0.38) Ref. [16] . (3.11) 

Only agrees in the two approaches. We have checked that scalar and tensor 
exchange, taken with the proper sign, generate values for a\ 2 that are not in 
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disagreement with Eq. (3.11) - as has been foreseen in the comments made in 
Ref. [16] concerning these two approaches. It would be very useful to recal- 
culate these couplings, by minimizing the amount of information used which 
goes beyond what is known from QCD, e.g., along the lines outlined in [17]. 
Finally, we note that in Ref. [18], C34 has been determined from a chiral sum 
rule. 

We now shortly discuss the uncertainties that we shall attach in the following 
to these couplings. In the case of 6^, we shall use 



6'' = 0.4 ±0.4. (3.12) 

As far as the polarizabilities are concerned, {cti—Pi).^± is independent of and 
determined precisely by the chiral expansion to two loops, once is fixed. We 
will then simply display this quantity as a function of al - the result turns out 
to be rather independent of its exact value, see Subsection 6.2 for a detailed 
discussion, and the uncertainty to be attached to it does not, therefore, matter 
here. On the other hand, those polarizabilities which depend on cannot be 
determined precisely from a calculation to two loops for reasons explained in 
Subsection 6.2 - we do not, therefore, worry here about the precise value and 
uncertainty for 

To complete this discussion, we note that we shall use = 92.4 MeV [19] (see 
[20] for a recent update of this value), and = 139.57 MeV in numerical 
calculations. 



4 Evaluation of the diagrams 



The lowest-order contributions to the scattering amplitude are described by 
tree- and one-loop diagrams. These contributions were calculated in [21]. The 
two- loop diagrams are displayed in the Figs. 1, 3 and 4. The two- loop diagrams 
in Fig. 1 may be generated according to the scheme indicated in Fig. 2, where 
the filled in blob denotes the rf-dimensional elastic TTTT-scattering amplitude at 
one-loop accuracy, with two pions off-shell. 

The diagrams shown in Fig. 3 may be reduced to tree-diagrams by using Ward 
identities [22]. They sum up to the expression 



2Z^g'^''-{{2p,-q,r{2p2-q2) 



1 



Ml-t 



Z^R{t) 



+ crossed > , (4.1) 



where is the pion renormalization constant. The function R{t) starts at 
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order and can be obtained from the full pion propagator [22]. 

Two further diagrams are displayed in Fig. 4. The first one - called "acnode" 
in the literature - may again be evaluated by use of a dispersion relation, see 
[5]. The second one is trivial to evaluate, because it is a product of one- loop 
diagrams. The remaining diagrams at order are shown in Fig. 5. 

The evaluation of the diagrams was done in the manner described in [5,23] 
and invoking FORM [24]. In particular, we have verified that the counterterms 
from the Lagrangian [12] remove all ultraviolet divergences, which is a very 
non-trivial check on our calculation. Furthermore, we have checked that the 
(ultra-violet finite) amphtude so obtained is scale independent. 




(1) (2) (3) (4) 




(5) (6) (7) (8) 




(9) (10) (11) (12) 




(13) (14) (15) (16) 




(17) (18) (19) 



Fig. 1. A set of two- loop diagrams generated by £2 and one-loop diagrams generated 
by £4- 
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Fig. 2. Construction scheme for the diagrams in Fig. 1. 
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Pi 



+ 



92 



P2 



Fig. 3. A class of one-particle reducible diagrams. The filled in circles summarize 
self-energy and vertex corrections. 





(1) (2) 

Fig. 4. Acnode and butterfly diagrams. 
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(1) (2) (3) (4) 



Fig. 5. The remaining diagrams at order p^: one- loop graphs generated by £4, and 
counterterm contributions from 

5 The two-loop amplitudes 



We give the expressions for the amplitudes A and B using the same notation 
as in [4]. This results in 



= + + ^ {g.(^) + ill} + + + °(V\ (5.1) 

The unitary part f/^ contains s, t and u- channel cuts, and P4 is a linear 
polynomial in s. Explicitly we find, 



2887r2sFi 



s - AMI) H{s) + 4 [ sG{s) + 2Ml{G (s) - 3 J (s)) ] 



00 



+ 



A^{s,t,u) , 



s - AMI) H{s) + 4 [ sG{s) + 2Ml{G (s) - 3 J {s)) ] 



00 



(5.2) 



with 



c{s, = ^ {ii - 1) iiQs' - hmls + UMt) 

+ {I2 - ^) (8^2 - 24M2s + 32M^) 
-12M^r3 + 12M>/4 - 12M> + 12M^} , 
4 = ^(3cos2^-l). (5.3) 

The loop functions J etc. are displayed in Appendix B, and Gt,{s) in Eq. (5.1) 
stands for G{s) evaluated with the physical mass. The term proportional to 
d^Q in Ua contributes to Z)-waves only. For see below. The polynomial part 
is 
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1 



(167r2F2)2 

„ 1 



( - - - 337 
ai = ai + -<!4/2 + / f -10 + 18 ^2 - 12 /a + — 

-^ri-5r2 + i2r4rA + 4| , 

If /I- 3- - 127 

5 - 5- - 21 
-Y^^i-4^2 + 3/A + y 



(5.4) 



The result for B reads 



with the unitary part 



1927r2sF4 1 3 
For the polynomial part we find 



{i (fi - i) + (^2 - ^) I + A.(s, (, u) . (5.6) 



(5.7) 

The integrals Ayi_B(s, it) contain contributions from the two-loop box, ver- 
tex and acnode graphs and also from the reducible diagrams. The explicit 
expressions for these quantities are given in Appendices C and D ^ . 

As an apphcation of the above, we plot the total cross section in Fig. 6, using 
the LECs from Eqs. (3.8) - (3.10). The data are taken from [25]. It is seen 
that the two-loop corrections are tiny in this kinematical region. 

In order not to interrupt the argument, a detailed comparison of our result with 
the previous calculation performed by Burgi [4] is relegated to Appendix A. 



(167r2F2)2 ' 

L z,r 1 f,2 , 3- 53\ 1 - 

^ = ^ -TsV +H2^^+2^^-24 -12^^- 



The corresponding Fortran codes are available upon request from the authors. 
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Fig. 6. The 77 — > tt'^tt cross section a{s; \ cos9\ < Z = 0.6) as a function of ^/s. 
The experimental data are taken from [25]. 

6 Pion polarizabilities: dipole and quadrupole 

The dipole and quadrupole polarizabilities are defined [26,27] through the ex- 
pansion of the helicity amplitudes at fixed t = M^, 



— t = M^) = (ai ± (3,U + tit(«2 ± P2U + 0{s') , (6.1) 



where denote the helicity amplitudes with Born-term subtracted. 
Because we have at our disposal the helicity amplitudes at two-loop order, we 
can work out the polarizabilities to the same accuracy. It turns out that all 
relevant integrals can be performed in closed form. We discuss the results in 
the remaining part of this Section. 

6.1 Chiral expansion 

Using the same notation as in [4], we find for the dipole polarizabilities 



TT 



12 




(6.2) 



12 



where 



2 - 

Cl+ = 0, Cl_ = g 



o,r 4 r, / 1 - 3 - \ 53 , 1 ^ 3 ^ 91 ^1 
+^ + A.}. (6.3, 

with 



8105 135 2 A 41 53 2 

7^^ A_ = TT^ . (6.4) 

576 64 ' 432 64 ^ ^ 



For the quadrupole polarizabilities, we obtain 



with 

C2+ = , C2- = 2 , 

2062 10817 2 8 - 8 - 

^^+ = "^ + 1440"^ +45^^+15^^' 

- S - - 218- 238 

d2- ^12al-2Ab'-l (10 + 4.h)-—l3+4.h-4.lA-— h - -rrh 

-55-™.^ (6.6) 
45 1920 ^ ' 

Wc end this subsection by evaluating the polarizabihties, using the above 
expressions and the central values for the LECs in Eqs. (3.8)-(3.10). The results 
are shown in Table 1 ^ . The numbers in brackets correspond to the order 
LECs in Eq. (3.11). The uncertainties are discussed in the following subsection. 



Dipole (quadrupole) polarizabilities are given in units of 10 ^ fm^ (10 ^ fm^). 
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Table 1 

The dipole and quadrupole polarizabilities. The numbers in brackets correspond to 
the order LECs in Eq. (3.11). 





to one loop 


to two-loops 


(ai - 


6.0 


5.7 [5.5] 







0.16 [0.16] 


{0L2 - (32) n+ 


11.9 


16.2 [21.6] 


{a2 + 02)^+ 





-0.001 [-0.001] 



6.2 Estimating the uncertainties 

The uncertainty in the prediction for the polarizabihty has two sources. First, 
the low-energy constants are not known precisely. Second, we are dealing here 
with an expansion in powers of the momenta and of the quark masses. We 
carried out this expansion up to and including terms of order p^. Higher order 
terms will influence the result - by which amount? 

We start the discussion by considering the uncertainties in the LECs. We shall 
use the order p'^ LECs displayed in Eqs. (3.8) and (3.9). The LECs at order p^ 
are not well known, see the discussion in Section 3. In Table 2, we display 
the contributions from the LECs at order p^ to the four polarizabilities, using 
the values from Eq. (3.10). The ones corresponding to Ref. [16] - displayed in 
Eq. (3.11) - are given in square brackets. The only significant difference occurs 
in the value of the difference of the quadrupole polarizabilities {a2 — (^2)-k^- 



Table 2 

The contribution of the LECs at order to the polarizabilities, according to 
Eq. (3.10). The numbers in brackets correspond to Eq. (3.11). 





a\ 


T 

«2 




total 


(ai - Pi)t,± 


-0.14 [-0.37] 





0.14 [0.13] 


[-0.24] 


(qi +Pi)^± 








0.14 [0.13] 


0.14 [0.13] 


("2 - P2)n± 





0.72 [6.09] 


-0.83 [-0.78] 


-0.10 [5.31] 


(a2 + /?2)7r± 















Resonance saturation generates a second inherent uncertainty - should one 
saturate at /i = 500 MeV or at // = 1 GeV? In Table 3 we display the 
polarizabilities for saturation at these two scales, using Eq. (3.10). It is seen 
that the induced change (which is independent of the values of the LECs at 
order p^) is quite substantial for (0^2 — f^2)-ir±- This is related to the fact that 
this quantity contains a rather substantial chiral logarithm, see below. 
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Table 3 

Resonance saturation induces a scale dependence in the amplitudes. Displayed are 
the values of the polarizabilities in case that saturation is assumed at ;u = 500 MeV 
or at /X = 1 GeV. 





II = 500 MeV 


fi = 1 GeV 


("1 - Pl)n± 


6.1 


5.5 


(ai +Pl)n± 


0.20 


0.13 


("2 - f32)^± 


14.6 


17.2 


(a2 + /32)7r± 


-0.001 


-0.001 



We now discuss the second source of uncertainties, the truncation of the chiral 
expansion itself. It is clear that, to have an idea of higher order terms, one 
needs at least the first two terms in the expansion. This makes it already 
clear that it is difficult to make rehable predictions for the polarizabilities 
connected with the hclicity flip amplitude, from which we have determined 
here the leading order contribution only. So, let us concentrate first on the 
helicity non-flip case 

The Born-term subtracted helicity amplitude does not have branch points 
at the Compton threshold. This is why it can be expanded there into an or- 
dinary Taylor series e.g. in the variable s and u = (t — u). One then expects 
that the amplitude at the Compton threshold is less affected by chiral loga- 
rithms than its counterparts at the threshold for 77 — > tttt, where unitarity 
cusps do occur. This is illustrated in Fig. 7, where we display the quantity 
10^M^iJ_i__i_(s, t = -u) as a function of s at t = u. Above the threshold s = 4M^, 
the modulus is shown. The solid (dashed) line is the expression to two loops 
(to one loop). It is clearly seen that the corrections at the Compton threshold 
are much smaller than the ones at the threshold for 77 — > tttt. To identify 
the chiral logarithms, we note that according to Eq. (3.7), the quantities k 
diverge in the chiral limit, li = —I + A[, where A[ is quark mass independent. 
We decompose the in the amplitude in this manner and display the 

contributions from the chiral logarithms / in Table 4. The numbers illustrate 
that, indeed, chiral logarithms generate a smaller contribution at the Compton 
threshold than at threshold for pion pair production. 

As for the quadrupole polarizabilities (0:2 — (32)tt±, also connected with the 
non-flip amplitude, it is seen from Table 1 that there is a substantial two- 
loop correction to the one- loop result. This can be again understood from 
the behavior of As its value is considerably increased at the threshold 

77 — s> TTTT, its slope at the Compton threshold receives a substantial correction 
as well, in order to make up that change. Indeed, the chiral expansion generates 
chiral logarithms that are responsible for the major part of the increase. If we 
decompose (0:2 — (32)-!r^ in a rnanner analogous the amplitude above, we 
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Fig. 7. The helicity non-flip amplitude in units of as a function of s, at 
t = u, with Born term subtracted. For s < 4M^ the quantity shown is 10^M^ff++, 
and for s > 4M^ we display 10^M^|ff++|. The solid (dashed) line is the expression 
to two loops (to one loop). The Compton threshold in jTr^ — > 77r^ and the threshold 
in 77 — > TT~^TT~ are denoted by the encircled letters C and A, respectively. It is clearly 
seen that two-loop corrections are suppressed at the Compton threshold. 



Table 4 

The amplitude W'^ M^H^^{s,t = u) at the Compton threshold s = and at 
s = 4M^. The contribution from chiral logarithms, listed in the fourth column, 
is included in the two-loop result quoted in column three. The normalization is 
Nh = lO^M^. The LECs at order are the ones from Eq. (3.10). 



NHH++{s,t = u) 


to 1 loop 


to 2 loops 


chiral logarithms 


s = 


2.89 


2.77 


-0.35 


s = ml 


7.13 


8.80 


1.28 



find that chiral logarithms contribute with ~ 4.5 x 10 ^ fm^ at two-loop order. 
These logarithms are, of course, independent of the LECs at order p^. 

Finally, we have checked whether there are potentially large contribution to 
at order . Using the same procedure as in [5], we found that all contri- 
butions from resonance exchange with masses below 1 GeV have a negligible 
effect - we do not quote the corresponding numbers here. 
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We now turn to the helicity flip amplitude which starts out at order p^: 
we have determined here only its leading order term in the chiral expansion. 
We checked whether there are potentially large contribution to at or- 
der p^, as is the case in 77 — > 7r°7r° [5]. Whereas, there is substantial contribu- 
tion from a;-exchange in the neutral case, this resonance does not contribute 
here, and the remaining contributions from p-exchange are very small, ex- 
cept for the contribution to the slope, parametrized by {0.2 + 02)it±, which 
is affected by —0.04 x 10"*^ fm^. On the other hand, there are also contribu- 
tions from one- loop graphs at order p^, where each vertex is generated by an 
anomalous contribution 771 — > tttt from the Wess-Zumino-Witten Lagrangian 
at order p^. We see no reason why these should be small compared to the 
leading term at order p^. We, therefore, do not consider the chiral prediction 
for these quantities particularly reliable. 



6.3 Values of the polarizabilities 

We have now all ingredients to provide a value for the polarizabilities in the 
chiral expansion, and start the discussion with (ai — Pi)jj.±. We add the uncer- 
tainties from the couplings at order p^, from and from the scale dependence 
introduced by the resonance scheme in quadrature and obtain 



for the so generated uncertainty. We display in Fig. 8 the result as a function 
of the rather inaccurately known coupling a^. We indicate with a filled square 
[filled circle] the two-loop result, evaluated with the LECs given in Eq.(3.10) 
[3.11]. The width of the band is twice the uncertainty displayed in Eq. (6.7). 
Let us note that the two-loop prediction differs only slightly from the one-loop 
calculation, see Table 1. This again shows that the value for the dipole polariz- 
ability is rather reliable - there is no sign of any large, uncontrolled correction 
to the two-loop result. We use the maximum deviation 1.0 from the central 
value 5.7 as the final theoretical uncertainty for the dipole polarizability, and 
obtain 



In the same figure, we also show the result of the recent measurement at 
MAMI [8] of this quantity. It is seen that the chiral expansion is in confiict 
with that measurement, independent of any reasonable value chosen for a^. 
The discrepancy with the recent dispersive analysis by Fil'kov and Kashevarov 
[10], displayed in Table 5, is even more pronounced. 



A = 0.80 X 10"^ fm^ 



(6.7) 
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Next, we discuss the quadrupole case (0:2 — ^2)^^- Here, the situation dif- 
fers. First, as we discussed above, chiral logarithms at order generate 
now a rather large contribution. As a result of this, the scale dependence 
of the saturation scheme is pronounced as well, see Table 3. Second, the 
LEG ttg enters with weight 12. As a result of this, it now matters which 
value of al is used, even though its contribution is suppressed by the fac- 
tor M2(167r2F2)-i ~ 0.014. Although the effect of using the [16] value for 
goes in the right direction and brings the result in closer agreement with the 
analysis of Fil'kov and Kashevarov [10], 

(a2 - P2U = (25.0t[J|) X 10-^ fm^ , (6.9) 

we cannot take the outcome as a reliable two-loop prediction of ChPT: com- 
pared to the one-loop result, the two-loop contribution generates nearly a 100 
percent contribution in this case, see Table 1. 

As we discussed in the previous subsection, the situation is no better for 
the polarizabilities (ai^2 + ^1,2)^^ associated with the helicity flip amplitude: 
considerably more work is required to put the chiral prediction on a firm basis. 

16 1 ' \ ' I I I I I I I 




A I I I I I I I I I I I 

-10 -8 -6 -4 -2 

r 

^1 

Fig. 8. The polarizabilities (ai — Pi)-,^^ at two-loop order, as a function of a^. The 
filled in square [filled in circle] denotes the value obtained by using for a[ the value 
from Eq. (3.10) [(3.11)]. The width of the ChPT band is calculated in the manner 
described in the text. "MAMI data" refers to the recent MAMI experiment [8]. We 
have added the errors quoted there in quadrature. 
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Table 5 

Experimental information on (qi — /3i)^±, in units of 10"*^ fm"^. We indicate the 
reaction and the data used. In [28] and [29], ai was determined, using as a constraint 
ai = — /?!. To obtain {ai — (3i)^±, we multiplied the results by a factor of 2. 



Experiments 


(ai - 


7P 77r+n Mainz (2005) [8] 


11.6 ± 1.5stat ± 3.0syst ± O.Smod 


L. Fil'kov, V. Kashevarov (2005) [10] 




77 ^ TT+TT- MARK II [25] , 




TPC/27 [30] , CELLO [31] , 




VENUS [32] , ALEPH [33] , BELLE [34] 




A. Kaloshin, V. Scrcbryakov (1994) [35] 


5.25 ± 0.95 


77 TT+TT- MARK II [25] 




Crystal Ball Coll. [36] 




J.F. Donoghue, B. Holstein (1993) [28] 


5.4 


77 TT+TT- MARK II [25] 




D. Babusci et al. (1992) [29] 




77^ TT+TT- PLUTO [37] 


38.2 ± 9.6 ± 11.4 


DM 1 [38] 


34.4 ±9.2 


DM 2 [39] 


52.6 ± 14.8 


MARK II [25] 


4.4 ±3.2 


7P ^TT+n Lcbcdcv Inst. (1984) [40] 


40 ±24 


TT'Z 7TT-Z Serpukhov (1983) [41] 


15.6 ±6.4stat± 4.4 syst 



Orders of magnitudes can be read off from Table 1. 



In Table 5, we collect experimental information on the dipole polarizabilities of 
the charged pions. In view of the large scatter of the results, which are partly 
inconsistent among each other, it is interesting to note that the COMPASS 
collaboration at CERN is presently remeasuring pion and kaon polarizabilities 
with high statistics, using the Primakoff effect [42]. It is also worth mentioning 
that a reanalysis of the MAMI data [8] , using a chiral invariant framework, is 
underway [43]. 



7 Summciry and outlook 

1. We have recalculated the two- loop expression for the amplitude 77 
TT+TT- in the framework of chiral perturbation theory. We have made use 
of the techniques developed in Ref. [5,23], and of the effective Lagrangian 
Cq constructed in [11,12]. 

2. The two Lorentz invariant amplitudes A and B are presented as a sum 
over multiple Feynman parameter integrals, whose numerical evaluation 
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poses no difficulty. This is in contrast to Rcf. [4], where part of the amph- 
tudes, denoted by A 4^3, could be published in numerical form only. Fur- 
ther, the method has allowed us to evaluate the dipolc and quadrupole po- 
larizabilities in closed form. As far as we are aware, the quadrupole polar- 
izabihties for charged pions have never been calculated in ChPT before. 

3. Om result agrees with the earlier calculation [4] up to the remainder 
Aa,b- The induced changes in the cross section and in the dipole polar- 
izabilities are far below the uncertainties generated by the (not precisely 
known) values of the low-energy constants. For the cross section below 
500 MeV, the change is less than 1 percent. For the polarizabilities, the 
change is given in columns two and three of Table A.l. 

4. We have investigated the uncertainties in the polarizabilities due to higher 
order corrections, and due to the uncertainties in the LECs at order p^. 
According to our analysis, the two-loop result Eq. (6.8) for the dipole 
polarizability (cti — /9i)^± is particularly reliable. It is in conflict with the 
recent experimental result obtained at MAMI [8], or with the dispersive 
analysis performed in [10]. 
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A Compcirison with the previous calculation 

In this appendix, we compare the amplitudes A, B and the dipole polariz- 
abilities with the previous calculation presented in Ref . [4] . In that reference, 
the amplitudes were evaluated with a different technique. Furthermore, the 
Lagrangian Cq was not available in those days, and an important ingredient 
to check the final result was, therefore, missing. We can make the following 
observations. 

(1) The amplitudes A and B consist of a part with explicit analytic expres- 
sions, and additional terms A^^^j that are displayed in Appendices C and 
D of the present work in the form of integrals over Feynman parameters. 
Aa,b were given only in numerical form in [4]. 
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(2) We find that our explicit analytic expressions agree with the ones in [4]. 
To compare A^^, we made two checks. First, we evaluated the cross 
section for the reaction 77 tt+tt" below a centre-of-mass energy of 500 
MeV, using the same values for the LECs as in [4] , 

{h, k, k, Ia) = (-1.7, 6.1, 2.9, 4.3, 2.7) , 

{a[, al, }f) = (-3.3, 0.75, 0.45) . (A.l) 

Our result agrees within a fraction of a percent with the two-loop result 
displayed in Fig. 9 of Ref. [4]. 

(3) Second, we evaluated the dipole polarizabilities by using the values of 
the LECs in Eq. (A.l). The result is displayed in the third column of 
Table A.l. In the second column, we display the values obtained in [4]. 
It is seen that the two results agree very well. The small difference is 
entirely due to the slightly different values of the quantities A^ ^- 

Table A.l 

Comparison of the polarizabilities with the previous calculation by Burgi [4]. The 
second column displays Burgi's result, the third one our evaluation, using Eq. (A.l). 
The fourth column displays the polarizabilities evaluated with the LECs used in the 
present work. 





Burgi [4] 


Present work 


Present work 






LECs from Eq. (A.l) 


LECs from Eqs. (3.8) - (3.10) 


(«! - /3l)^+ 


4.42 


4.39 


5.72 


(ai +/3i)^+ 


0.31 


0.28 


0.16 



The following comments concerning the central values of the polarizabilities 

are in order. We display in the last column in Table A.l the values of po- 
larizabilities obtained by using the LECs from Eqs. (3.8) - (3.10). It is seen 
that the combination (ai — I3i)t^± differs considerably from the one obtained 
with Eq. (A.l). To identify the source of this difference, we first note that the 
polarizabilities depend linearly on the LECs, except the term l/^l^ in as 
can be seen from the Eqs. (6.2-6.4). We then expand {ai — (3i)t^± in the LECs 
around their central values, drop quadratic terms and obtain at = Mp the 
decomposition 

(ai -/3i)^± =5.719 

+ 0.215 ■ {h + 0.4) - 0.177 ■ {I2 - 4.3) - 0.014 ■ ([3 - 2.9) 

+ 0.172 • ([4 - 4.4) + 2.44 ■ (/a - 3) + 0.04 • {a\ + 86' ) . (A.2) 

Inserting in this expansion e.g. the values in Eq. (A.l) generates the result in 
column 3 of Table A.l. It turns out that all changes in the order LECs have 
conspired to modify Burgi's result for (cti — /9i)^± towards positive values. The 
contribution from the LECs at order is negligible. 
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In summary, we conclude that our calculation nicely confirms Burgi's result 
[4], provided that the same set of LECs is used. 



B One-loop functions 



In order to simplify the expressions, we set the pion mass equal to one in the 
following Appendices, 



(B.l) 



1. The loop- integral G{s) is 



(B.2) 



G is analytic in the complex s - plane, cut along the positive real axis for 
Re s > 4. At small s. 



(n!) 



1^2 



-j^ oo 

^^^^ " 16^ + l)(2n + 1)! 



(B.3) 



The absorptive part is 



ImG(.) = J-ln|i±^|, 



s > 4, 



a = Jl-A/s, 



(B.4) 



and 



- 16t:^G{s) = < 



1 -iarctg^(^)2; 



4 < s, 
< s < 4, 
s < 0. 



(B.5) 



In the text we also need 



G{s)=G{s)-sG'{0). 



(B.6) 
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2. The loop- integral J{s) is 



J{s) 



167r2 



dxln{l — sx{l — x)) . 



(B.7) 



J is analytic in the complex s - plane, cut along the positive real axis for 
Re s > 4. At small s, 



1 °° CnM^ 



167r2 n{2n+l)\ ' 



(B.8) 



The absorptive part is 
a 



Im J(s) 



IGtt ' 



s > 4 



(B.9) 



Explicitly, 



167rV(s) = < 



2-2(^)^arctg(^)^; 



In the text we also need 



4 < s, 
< s < 4, 
s < 0. 



(B.IO) 



J (s) = J(s) - s J'(0) . 
3. The loop- function H is defined in terms of G and J, 
= (s-10)J(s) + 6G'(s). 



(B.ll) 



(B.12) 



C The quantities and 



Here we display the expressions for the quantities Aa{b)- 



AA{s,t,u) = 



1 f /35947 TT^X 1130291 
(47rF^)4 { \ 8100 ^ TJ ~ 6350400 



(47rF^)4 288 
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A / X 1 [ 492197 /81101 TT^ 1 ^ ^ 1^, x\l 
Afi s, t, u) = <^ + -S(t) + -S(u) - 



Here, the function S{t) is defined by 



SW^jML. (C.3, 



CXJ 1 

= - J da(5{a) j dx [Qt + (1 - I2x + ISx^jt^} 

4 

L ^2(a^,t) ^ (t-l)a;(l-a;) | ^ 
\ Z2(a:^, 1) ^2(3:;, 1) J ' 



/?((7) = Vl-4/(7, 

2;2(x, = + (7(1 - x) - - l)x(l - x) . 

The loop function hpit) stems from the sunset diagram. The functions Fi are 
stemming from the acnode, vertex and box diagrams and are defined by 



r r \ P 6 p 

Fr"= da3 dPxl ^^"^"^^^ ^^^pw 

^ ;^ yy ■ ^acn(l) ^=2 ^acn(i) 

oo 1 

4 

oo , ^1 2 

^box ^'^EHrbU^£l + ^Sx_^£L} , ^ = (C.4) 

/I n n— 1 



and 



r)(") _ - _|. 



1 1 



box;t box;u 

Here Pj are polynomials m. s^v — t — u and in Xj. Their explicit expressions 
are given in Appendix D. The arguments of the functions are defined by 



z^x.(i) =y- ttiX^il -X3), (i = 1, . . . , 6) , y = xl + a{l- X3) , 
Oi =02 = X1X2S + xi{t - 1) + X2{u - 1) , 
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a^ — t — 1 , — xi{t — 1) , 
a^^u-l, aQ = X2(u-l), 

Zyer = Ct(1 - X3) + xly2 , y2 ^ 1 - SX2{1 - X2) , 

^box;t — Bt — AtXi , 

A = X2X3 

Bt = At -\- ^ver ) 
^box;u ^box;t 1 1— >u ■ 



S{1 - X2)X3 + {I - t){l - X3) 



(C.5) 



The acnode integrals are easy to evaluate numerically in the physical region for 
the reaction 77 Tin, because branch points occur at t = 4, -u = 4 only. On 
the other hand, the vertex and box integrals contain branch points at s = 4. 
In order to evaluate these integrals at s > 4, we invoke dispersion relations in 
the manner described in [23]. 



D The polynomials Pa and Pb 



Here, we display the polynomials Pa{b) that occur in the expressions A.a{b) in 
Appendix C. We use the abbreviations 



— Xi ~\~ X2 — ^X\X2 ; X^ — X-\_ — X'2 , 

2^123 = (1 + 2^3 - 2X2X3) (1 - 2:3 + 2X1X2X3) . 



(D.l) 



D.l The polynomials Pa 



A; axMi(l) 



-PA;acn(2) 



+ 
+ 



144X3(1 — X3)(SX_(- — i^X_), 

12s^(l — X3)^|x+(5x+ + 6x^ ~ 1) ~ 5x^X3 — 6x^X3 
x^(9 — 6x+) — x+(2 — 7x+) x| 



6xi(x+ — 1) + (3 — 7x+)x+jx3| 
+12si/(l - X3)^x_{l - 6x^ - 14x+ + IOX+X3 



2x^ - 



2(6x? 



+6(1 

+ (6x^ + 8x+-3)x^} 
-12z/2(l - X3)2{6x+(1 - 2x+)x: 



l)x 



— X 



9 - 5x3 - I2X+X3 - (6 - 18x+)x3 + (5 - 12x+)x3] } 
+24i/(l - X3)x_{l + 2x+ - (25 + 2x+)x3 + 3(7 - 8x+)x3 
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+{100x+ - Al)xl + (51 - n8x+)xl + (42x+ - 19)x^} 
-12s(l - X3)[2x+{1 + 2x+) + 2(2 - 12xl_ - 17x+ - 2x^)^3 
+2(13 - 6x^_ - nx+)xl - 2(41 - 18x1 - S9x+ + 8x^)x| 
+ (82 - 36xt - 194:X+ + 28a;^)a;^ 
-2(15 - 6a;?. - 35a;+ + 6xl)xl} 

+48x3(1 - a;3){ - 2 - 4x+ + (35 - 50x+)a;3 

+71(2x+ - l)xl + (61 - mx+)xl + 2(21x+ - 10)xl} , 

-PA;acn(3)= 3{s - {1 ~ Xsf [1 + X3) 

+6(s - z/)(l - X3)^(3 + 5x3 + 2x1 - 2x1) 
+ 12x3(1 - X3)(2 - X3)(3 + 3x3 + xl- xl) , 

i"A;acn(4) = 24(s - ufxl{l - 2Xi)(l - X3)'(2 + X3) 

+48(s - t/)xi(l - 2xi)(l - X3)'(2 + 4x3 - 3x^) 
+96(1 - 2xi)x3(l - X3)(2 - X3)(2 + 2x3 - x^) , 

Pa; acn(5) — Pa; acn(3) | t<->u ) 

Pa ;acn(6) — -PA;acn(4) |t-<-»u, xi*-+X2 ; 



PA;ver = -8s'(l ^ 2X2)x2x^{ - 54 + 72(3 + 4X2)X3 

+ (330 - 1134x2 - 811x^)x^ - 108(5 - 15x2 - 4x^)x^ 
-45x2(18 - 75x2 + 32(3 - Xajx^jx^} 
+32s(l - 2x2)x^x^{33 + 8(27 - 7x2)x3 
-3(222 - 30x2 + 55x2)x3 - 540(1 - 3x2 + xl)xl 
+1215(1 -X2)'x^}, 



A;box;+ 



4s' 



|6xi 



9 - 2(23 + 8x2)x3 - (67 - 405x2 + 31x^)x^ 



+(70 + 39x2 - 808x^ + 20x 
+9(6 - 61x2 + 53x^ + 60x^)x^ - 81x2(3 - 10x2 + 8x^)x^ 
+3x^x2X3 — 92 — 171x3 + 592x2X3 

+ (231 + 4(194 - 339x2)x2)x^ 
+9(1 - 2x2) (71 - 22x2 - 60x^)x^ 
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+27(1 - 2x2)^(11 - 16x2)xl\ - 2xlxlxl[2A5 + (1 - 2x2)x3 
x(470 + 27(1 - 2x2)2:3(15 + 8(1 - 2x2)0:3))" 
+6a:3(-19 + 24x3 + 35x^ - 36x^ 
+X2X3(-25 + 40x3 + 27(5 - 6x3)x3) 
+X2(10 + 44x3 - 199x^ + 9x^(8 + 9x3)))} 

+ 12u^xlxj{ - 20 - 16xi(2 + 7xi) 

+ [48 + 50x2 + 310x1 + lOlx^ + 30xi(3 + 15xi - 5x?)x2 
+ f - 78 + 2xi + 411xi - 20(9 + 63xi + 57xi + 13x?)x2 



+120xi(l + xi + 5xi)x2 



x^-3 



24 + 3x1 (94 + 73xi) 



-2(69 + 351x1 + 381x? + 19x?)x2 
+4xi(24 + 117x1 + 34x?)x^]x^ 

+9(1 - 2x2) [I8 + 54x1 - 36xiX2 + x^(33 - 2(9 + 8xi)x2)]x^} 
-485x2X3! - 6 + [18 + 20x2 - 31XiX2]x3 
- [29 + 49x2 + 4x2 + 4x1X2(9 - (22 - 8xi)x2) 



Xn 



+ 



6 + 202x2 



54x2 + 164xiX2 



+ 114xiX2 



4x1X2 



xfx^(149 - 58x2) 



+ 



65 - 333x2 - 40x^ + 189xiX2 - 4xi(137 + 75x2)x^ 



+5x1X2(23 + 12x2(5 + X2)) 



X q 



-9 



4 - 46x^ - xfx^(l - 2x2) (41 + 8x2) 



76x^) 



+XiX2(25 + 18x2 
+81x2(1 - 2x2) 2 - 3Xi + Xi(2 + Xi)x2 - 2xiX2]x3} 
+96x2X^{ - 21 + [36 + 65(1 - Xi)x2 
+ [23 - 8(23 - xi)x2 + 176x1X2] X3 



2:3 



6 - (83 + 37xi)x2 + 120x1X2 X 



4 



29 + 10(17 - 7xi)x2 - lOOxix^ 
8 - (31 + xi)x2 + 32xix^]x^ + 81(1 -2c2)(l - 2xiX2)x^} , 



+ 



-8si/X2.X3|27xi — 
72 - 129xi + 3(6^ 



57 + 6x1(13 + (26 + 51xi)x2) 
+ xi(456 + 25xi))x2 



2^3 



+xi(282 + xi(1383 - 470xi))x^ 
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+ [105 + 330xi - 714x2 - 3xi(200 - 489xi)x2 
-2(75 + a;i(1557 + 912xi + 470x?))x^ 

-120xi(3 + 9x1 - 16x?)x^]x^ + 9[2x?x^(l - 2x2) (13 + 34x2) 
+3x^X2(1 - 2x2) (1 - 23x2 + 4x^) 



-6x1(1 - 2x2) (3 + 36x2 + 17x^) - 6(2 - 2x2 - 19x^) 
+27x2(1 - 2x2) [is + 18x1X2 - 33x? 
+2x^x2(33 - 9x2 + 8x1(1 - 2x2))]x^} 
+A8uxlxl[l8 - 29xi 

21 + 4(10x1 + 12x2 - 19(2 - xi)xiX2) 



X3 



+ 



42 + 106x2 + 252xi - 23x^x2(7 - lOxs 



-6x1X2(3 + 20x2) 



X3 



83 + 230x2 - 41xi 



+5x1X2(134 - 20x2 - 3xi(9 + 10x2)) 



-9 



12 - 54x2 + xi(l - 2x2) (53 - (24 + 41xi)x2) 



X Q 



+81(1 - 2x2) [2 + 3xi - xi(4 + xi)x2]x^} , 



= — 6s^ 1+Xi(l — 2x2) a;2^3a:;i23 

X 8 — 5x2 + 8x1 — xi(6 + 5xi)x2 — 6x3 

-6X1(1 - 2X2)X3 + 9X2(1 + Xi(l - 2X2))^X^ 
-6SZ/^(1 - Xi)^X^x|xi23 

X 8 - 6x3 + X2(l + xi(l - 2x2))(5 - 3x3(4 - 9x3)) 

+ 125^X2X3X12312(2 + X2 + 2Xi - (6 - Xi)XiX2) 



-12 



1+Xi(l-2X2)]X3+ [l3(2-X2) 



+2xi(l - X2)(13 - 16x2) - XiX2(13 - 32x2 + 32x3) 



Xn 



-12 



1 + Xi + X2 + X^(l - 2X2)^X2 - 4xiX2 



Xn 



+27x2[l+xi(l-2x2)] x^} 

-12i/^(l - xi)^x^X3Xi23{2 + 12x3 - 17x3 + 24x^ - 27x3} 

+24SX2X3X123 X |6(2 - X3)(l - X3)X3 

+2xix^ [23 - 24x3 + xl{2- 3x3) (2 - 9x3) 

+ (1 + xi)x2 [ - 23 + 24x3 - 4x2 + 3(8 - 9x3)xi] } 

+144x2(1 - X3) 2X3X123 5 + X3(2 + 3x3) , 
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p(2) 
■'A; box; 



6z/'^(l — a;i)^X2X3Xi23 — 6x3 + 9x1 
+6s^u(l — Xi)a:2X3Xi23 x |l6 + 16xi — 5x2 — 22xiX2 — 5x^X2 



2 + (1 + xi)x2 - 2xixl\ [1 + xi{l - 2x2) 



+27x2 



1 



(l-2a;2)]'a;^} 



— 24sz/(l — a;i)a;2a;3Xi23 x |2 — 6x3 + (13 — 6a:3)a;3 
-X2[l +xi{l- 2x2)] [6 - 2x3 + 9x^(2 - 3x3)] a;3} 
+72u(l — xi)(l — a;3)x2X3Xi23 5 - 7x3 + (1 - 9x3)x3] . 



D.2 The polynomials Pb 



B; acn(l) — 



-Pb; acn(2) ' 



-Pb; acn(4) 



144X3(1 - X3) 

12^/2 



x+ x. 

s 



+ 



s 

2A1/X- 



(1 — X3) (3 — 5X3 — X3)x_ + 6X3X+ 



(1 - X3)'(7 - 12x3 + 3x1 + 4^3 - ^4) 



48 



-— X3(l - X3)(14 + 19x3 - 43x1 + 29x^ - lO^a) 

+12s(l — X3)2|x_|_(l + 3x+) — 5x3X^ + GxgX^ 

— X3 9x1 — x+(2 + 3x+) + X3 6x1 — x+(3 + x+) | 

— 12i/x_(l — X3)2 1 + 6x+ — 10x3X4- + 6x^(1 + x+) 

+2x^(1 - 3x+) - x^(3 - 4x+)" 

-24(1 - xs)^ [7x+ - 2(1 + 6x+)x3 

-(1 - x+)xl(15 - 26x3) - (15 - 13x+)x^l , 



B; acn(3) 



3(s - z/) 



-(1-X3)'(1 + X3) 



-^^^^(1 - X3)'[3 + 5X3 + 2x1(1 - X3) 



12 

+— X3(l - X3)(2 - X3) [3(1 + X3) + xl(l - X3) 



24(s - 



■xl{l-XsY{2 + xs) 



29 



48(s - u) 



96 r 

X3{1 - X3){2 - X3) 2 + (2 - X3)X3 

s 



Pb; acn(5) — -PB;acn(3) | 



t<->u ) 



-PB;acn(6) -^B; acn(4) |t^u,xi<->X2 ) 



B; ver ■ 



85x2(1 - 2x2)xj[54: - 8(69 + 50x2)x3 



+15 
-180 



70 + X2(94 + 77x2 



X3 



3 + X2{4: + 23x2)]xl - A05x2{2 - 9x2)xt} 



-96x^(1 - 2x2)xl 25 - 144x3 + 500x^ - 780x^ + 405x^ 



p(i) 

B;box;+ ' 



12i/2 



-x^(l - X3)x^{20 + 16xi(7 + 9xi) 
s 

-98x3 - a:;i(574 + 573xi)x3 



+6 



288 



32 + xi(166 + 117xi)]x2 - 27[6 + xi(18 + llxi)]x^} 

X2(l - X3)^x^{7 - 17x3 + 14x2 _ -^7^3 ^ 27x^} 
18xi + [38 + 4(51xi - 5x2 + 45x1X2)1x3 



156 - 50x2 + 322x1 + 586xiX2 + 5x^X2(85 + 74x2) 
+5 [38 + 18x2 + 4xi + 278xiX2 + 76xiX2 



-3x^X2(7-4x2(17 + 7x2)) 



24 + 94x2 - 36xi + 130xiX2 + 424xiX2 



-3x^x2(1-2x2) (37 + 76x2) 



+27x. 



6 - xi(l - 2x2)(18 - llxi(l - 2x2))]x^} 



+48x2X3(6 - (34 + 25xiX2)x3 
+ [97 -X2(7- 16x1(7 + 4x2))" 
-2 [79 + X2(13 + xi(191 + 118x2)) 
+5 [25 + X2(45 + xi(101 + 110x2)) 
-3[l2 + X2(118 - xi(7 - 284x2))]x^ 
+81x2 [2 - 3xi(l - 2x2)] x^} , 



X Q 



X Q 



30 



—xl{l - X3)xlh4: + Uxi - 9{7 + 9xi)x3 
s 

+ (127 + 373xi)x2 - 12(19 + 46a;i)a;^ + 81(2 + 3a;i)a;^} 

+24uxlxl[ - 9xi + [l9 + 2a;i(41 + (38 + 81xi)x2)]x3 

- 78 + 34x2 + 189x1 

+xiX2(524 + 84x2 + 3xi(193 + 59x2))] Xg 

+5 [l9 + 38x2 + 34xi + xiX2(224 + 78x2 + 9xi(9 + 20x2)) 

-3 [2(6 + 53x2) + 18xi 

+XiX2(240 - lllxi + 248x2 + 447x1X2)1x3 



+27x2(6 + 18x1X2 - llxi(l - 2x2)] X3} , 



12i/2 



(l-Xi)^X^(l-X3)^X^ 



X 



2 + 12x3 - 17x^ + 24x^ - 27xi 



144 



X2(1-X3)^X^[5 + X3(2 + 3X3) 



+12sx^(l - X3)^x^{2(2 + X2 + xi(2 - (6 - Xi)x2)) 



-12 



1 + xi(l - 2x2)] X3 + [l3(2 - X2) + 26xi 



-XiX2(58 — 32x2 + xi(13 — 32(1 — X2)x2)) 



-12 
+27x2 



1 + Xi + X2 + xl{l - 2X2)^X2 - 4xiX2 



X., 



1 + Xi(l-2X2)] x^} 
-6s^x^(l + Xi(l - 2X2))(1 - X3)^X^ 

X {8 - 6x3 - X2(5 - 9x3) - XiX2(5 - 9(1 - 2x2)^X3) 
+2xi [4 - 3x3 - 18x^x^ - 3x2(1 - 2x3 - 3x^)] } 



-6z/^(l -Xi)^X^(l -X3)^X? 
xjS - 6X3 + X2[l + Xi(l - 2X2 
+6X2(1 - X3)'x^{24(2 - X3)(l - X3)X3 

-4(1 + xi)x2[23 - X3(24 - X3(2 - 3x3) (2 - 9x3)) 
+8x1X2 [23 - X3(24 - X3(2 - 3x3) (2 - 9x3))] } , 



5-3x3(4-9x3)]} 
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3(2) 
B;box:- 



(1 - xi fxl{l - X3)'^xl{5 - 6x3 + 9x1) 



2^6/ 



6z/^ 

s 

I (1 - xi)xl{l - xsfxl 5 - X3(7 - X3(l - 9x3)) 



+6sz/(l — Xi)x2(l — X3)^X3|l6(l + xi) — 5x2 — 22x1X2 



5x^X2 — 6 2 + (1 + Xi)x2 — 2x1X2 1 + Xi(l — 2X2 



X3 



+27xs 



1 + Xi(l-2X2)] x^} 
— 24z/(l — Xi)x2(l — X3)^X3|2 — X3 6 — X3(13 — 6X3 
+ (1 + Xi)x2(6 - X3(2 - 9(2 - 3X3)X3)) 

-2xix^(6 - X3(2 - 9x3(2 - 3x3)))] } . 
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